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Abstract The jets of microquasars with high-mass stellar companions are exposed 
to the dense matter field of the stellar wind as well as to the photon densities found 
in the surrounding medium. Photopion and proton-proton interactions could then 
lead to copious production of neutrinos. In this work, we analyze the hadronic mi- 
croquasar model, particularly in what concerns to the neutrino production. Limits 
to this kind of models using data from AMANDA-II are established. New con- 
straints are also imposed upon specific microquasar models based on photopion 
processes. These are very restrictive particularly for the case of SS433, a micro- 
quasar for which the presence of accelerated hadrons has been already inferred 
from iron X-ray line observations. 

Key words: X-rays: binaries — Stars: winds, outflows — gamma-rays: observa- 
tions — gamma-rays: theory 



1 INTRODUCTION 

The presence of relativistic hadrons in microquasar jets like those of SS433 has been already 
inferred from iron X-ray line observations (e.g. Migliari et al. 2002). A significant content of 
relativistic hadrons in microquasar jets could open the possibility for a hadronically-generated 
emission of high energy radiation. Here, we discuss a new mechanism for the generation of high- 
energy gamma-rays in microquasars that is based on hadronic interactions occurring outside the 
coronal region (Romero et al. 2004). In this model the gamma-ray and neutrino emission arises 
from the decay of neutral pions created in the inelastic collisions between relativistic protons 
ejected by the compact object and ions in the stellar wind. The only requisites for the model 
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to operate are a windy high-mass stellar companion and the presence of multi-TeV protons in 
the jet, both of which seem natural in microquasar environments. We pay particular attention 
to the possible neutrino signal of this kind of models, and impose constraints using the latest 
AMANDA-II data (Ahrens et al. 2004). 



2 THE HADRONIC MICROQUASAR MODEL 

2.1 The jet and the particle spectrum in the lab frame 

For simplicity, we shall not make any specific assumption about the magnetic field or other 
parameters in the jet, but rather model it as a beam of energetic particles (i.e., in the spirit 
of, for instance, Purniohamniad and Saniimi 2001, and Bednarek et al. 1990). The jet axis, 2, 
is assumed normal to the orbital radius a. We shall allow the jet to expand laterally, in such a 
way that its radius is given by R{z) = ^z*^, with e < 1 and zq < z < Zmax. For e = 1 we have a 
conical beam. The jet starts to expand at a height zq ~ a few hundred km above the black hole, 
outside the coronal region. The particle spectrum of the relativistic e ~ p flow is assumed to be 
a power law iV^ ^{E'^^ p) = K^, p E''^^^, vahd for E'^ p™'" < E'^ ^ < E'^^ p™'"'', in the jet frame. 
The corresponding particle flux will be Jg p{E'^^ p) = (c/47r)iVg p(i?e, p). Since the jet expands, 
the proton flux can be written as: 

4K) = ^^o(f)^"i?r^ (1) 

where n > (a value n = 2 corresponds to the conservation of the number of particles, see 
Ghisellini et al. 1985), and a prime refers to the jet frame. Note that these expressions are 
valid in the jet frame. Using relativistic invariants, it can be proven that the proton flux, in the 
observer (or lab) frame, becomes (e.g. Purmohammmad & Samimi 2001) 



Jp{Ep,0} = — J 



47r 



sin2 6' + r2 cos 61- 



0b Ep 



-,1/2 



(2) 



where F is the jet Lorentz factor, 9 is the angle subtended by the emerging photon direction 
(assumed to be similar to the initial proton direction) and the jet axis, /3b is the corresponding 
velocity in units of c. Note that only photons emitted with angles similar to that of the incli- 
nation angle of the jet will reach a distant observer, and thus 9 can be approximated by the jet 
inclination angle. 

In order to obtain Eq. ^ one has to consider conservation of particles. If n is the number 
of protons per unit energy per unit solid angle per unit volume, so that in the frame co- 
moving with the jet one has n'{E',n') = {A/ATT)E'-°'dV' , the equality n{E,n)dVdQdE = 
n'{E',n')dV'dn'dE' holds, what implies n{E,VL) = n' [E' ,n'){dn' / dVL){dE' / dE){dV' / dV). 
From the invariance of d^p/E, the invariance oipdQ.dE is also proved (Hayakawa 1969, p. 715), 
so that {dil/dn')(dE' /dE) — pjp' . Using relativistic Lorentz transformations, p'^ = F(p|| -|-/3i?), 
with pj| = pcos,9 and p'^ — P^ — psin6', what implies p'^ = p^(ii\Ti'9 -f- F2(cos6' + 

F and 



PE/VW^ 
E' = T{E 



m^)^). This defines p/p', which, together with the equalities dV'/dV 
- (3VE^ - w? cos 6*), yields to Eq. @. 



We will adopt the jet-disk coupling hypothesis proposed by Falcke & Biermann (1995) 
and applied with success to AGNs, i.e. the total jet power scales with the accreting rate as 
Qj — gjAfdiskC^, with gj = 10^^ — 10^^. The number density no' of particles flowing in the jet 
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at -Ro = R{zo) is then given by cttRqIIo' = QJuipC^, where nip is the proton rest mass. This 
impUes: 

Additionally, n'^ = ^^ZI N'p{E'p, zq) dEp . Then, if E"^^"" » E'f'', which is always the 
case, we have 

Ko = no'ia-l){E'^'^r-\ (4) 

which gives the constant in the power-law spectrum at zq- This completely defines the proton 
spectrum. 



2.2 Simple wind modelling 

The structure of the matter field in the wind will be determined essentially by the stellar mass 
loss rate and the continuity equation: M* = ATrr^p{r)v{r), where p is the density of the wind 
and V is its velocity. Hence, 

The radial dependence of the wind velocity is given by (Lamers & Cassinelli 1999): 

v{r)=v^{l-^fy , (6) 

where is the terminal wind velocity, is the stellar radius, and the parameter /? is 1 

for massive stars. Hence, using the fact that = + and assuming a gas dominated by 
protons, we get the particle density of the medium along the jet axis: 



n(z) = ^ 1 - , * . (7) 



Typical mass loss rates and terminal wind velocities for O stars are of the order of 10~^ A/q 
yr~-^ and 2500 km s~^, respectively (Lamers & Cassinelli 1999). This simple modelling for the 
wind was also used when analyzing the possible TeV emission from stellar systems (Romero & 
Torres 2003, Torres et al. 2004). 

It is important to note that wo arc considering not that the beam interacts with the wind 
in a face-on collision, but instead, that the wind diffuses into the jet from the side. The wind 
penetration into the jet outflow depends on the parameter w ~ vR{z)/D, where v is the 
corresponding velocity of wind/beam in the direction of diffusion, R{z) is the radius of the jet 
at a height z above the compact object, and D is the diffusion coefficient, w measures the ratio 
between the diffusive and the convective timescale of the particles. In the Bohm limit, with 
typical magnetic fields So ^ 1 — 10 G, < 1, and the wind matter penetrates the jet. 



2.3 Gamma- ray and neutrino emission 

The dominant 7r-producing channels in the hadronic interactions of the jet with the wind are 

(iJth ~ m^c^): 

p + p P + tt" ^ p + 2j (8) 

p + p ^ n + tt'^ ^ p + e~ + i>e + e'^ + t^e + t^n + J^fi (9) 
p + p — > p + 'JT~^+n~ 

^p + e+ + + +^1^1, + e-+Pe + '^J^n, (10) 
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where all symbols have their usual meaning in particle physics. Processes and (|1(JII lead 
to in situ pair creation and neutrino production. The differential gamma-ray emissivity from 
7r°-decays is: 



qy{Ey) = 4TT(7pp{Ep)- 



2Z 



(") 



JpiEj, 0)r]A- 



(11) 



Here, the parameter rjA takes into account the contribution from different nuclei in the wind 
and in the jet (for standard composition of cosmic rays and interstellar medium rj^ = 1.4 — 1.5, 
Dormer 1986). Jp{Ej) is the proton flux distribution evaluated at E = E^. The cross section 
(7pp{Ep) for inelastic p — p interactions at energy Ep « lOE^ can be represented above Ep « 10 

GeV by app{Ep) « 30 x [0.95 + 0.06 log(£;p/GeV)] mb. Finally, Z^^^^o is the so-called spectrum- 
weighted moment of the inclusive cross-section. Its value for different spectral indices a is given, 
for instance, in Table Al of Drury et al. (1994). Notice that is expressed in ph s~^ erg~^ 
when we adopt CGS units. 

The spectral gamma-ray intensity (photons per unit of time per unit of energy-band) is: 



I^{E^,e) = / n{r')q^{r')(f'r' 
Jv 



(12) 



where V is the interaction volume. The spectral energy distribution is L'!^ {E.y, 9) = E'i^I^{E^, 9) 
and using eqs. ©, HI), 0, dJ and lO, we get: 



Ij^O ^p^ttO a - 1 . ,mins„_i 

[E p ) X 
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r-"+i (E^- (3b J E^-mlc'^ cos i 
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dz. 



(13) 



This expression gives approximately the 7r°-decay gamma-ray luminosity for a windy micro- 
quasar at energies E^ > 1 GeV, in a given direction 9 with respect to the jet axis. The neutrino 
spectrum roughly satisfies (e.g., Dar & Laor 1997): dF^/dE^ ~ 0.7dFj/dE^. 



3 DISCUSSION ON NEUTRINO UPPER LIMITS 

Upper limits imposed by AMANDA II (Ahrens et al. 2003) in known microquasars are typically 
F < 4 X 10"''(i?/GeV)~2GeV^^cm~2s~i. This upper limit and the theoretical expectation for 
the neutrino production, given assumptions on the free parameters of the system, allows already 
to put some constraints on the most favorable hadronic models. For instance, the model in which 
the distance is D — 5 kpc or less, qj—10~'^ or higher, the stellar mass loss rate is 10~^ Mq yr~^ or 
higher, and the jet is inclined only 10 degrees or less, is already ruled out by neutrino astronomy. 
Note that this model maximizes the gamma-ray and neutrino emission, by enhancing the target 
matter density for proton interactions, by locating the microquasar at roughly half the distance 
to the Galactic Center, and by assuming a matter content in the jet with a significant dose of 
hadrons. In any case, this is a strong result in the sense that there is no need for any pointing 
of any telescope (the only assumption is that the microquasar is in the northern hemisphere, in 
the field of view of the AMANDA-II experiment). Note, however, that still there are plenty of 
models (different choices of system parameters) between the atmospheric neutrino background 
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Table 1 Comparison between upper limits and predictions of the photopion model 
for some microquasars. 



Microquasar 


Prediction 

(Distcfano ct al. 2002) 

erg cm~^ 


Upper limit 

(data from Ahrons et al. 2003) 
erg cm~^ 


SS433 


1.7E-9 


2.6E-10 


Cyg X-3 


4.0E-9 


1.3E-9 


Ci Cam 


2.2E-10 


2.9E-10 



and the typically imposed upper limits; these have improved observational expectations for 
ICECUBE (as an example, models with proton slope 2.2 are shown in Figure 1). 

Note that alternative hadronic microquasar models make use of photopion production of 
neutrinos (Levinson & Waxman 2001). The imdcrlying idea in this kind of scenarios is that 
neutrinos are the result of photopion processes with synchrotron photons and external target 
fields surrounding the microquasar jet. Distefano et al. (2002), have presented the expectations 
for the neutrino fluxes from known microquasars under the assumed validity of this model, and 
focused on their scrutiny with the forthcoming ICECUBE neutrino telescope. Apart from a 
technical point -ICECUBE search bin is most likely going to be 1 degree, not 0.3, what would 
enhance the background signiflcantly from Distcfano ct al. estimates, making detection more 
difficult- the strong claim in their study is that they predict in some cases even hundreds of 
events per year in a km-scale detector like ICECUBE. These would make of microquasars the 
most notable neutrino sources in the sky. For instance, in the case of SS433, the prediction 
is 252 muons events/yr above 1 TeV. But this implies around 20 or more events per year in 
AMANDA-II, what is ruled out already by data. Table 1 shows a comparison between the 
predictions of Distefano ct al. (2002) and the measured upper limits extracted from the work of 
Ahrens et al. (2004), assuming a spectrum, for some of the cases that appear to be already 
ruled out, or are on the verge of being ruled out, by AMANDA-II data. Note that the case of 
SS433 is especially significant: in this case it is known for sure that in the microquasar's jets 
there are protons and heavy nuclei accelerated to relativistic speeds. For SS433, the estimations 
of Distefano et al. are ruled out by about one order of magnitude. Further analysis of neutrino 
upper limits will be presented elsewhere. 



4 CONCLUDING REMARKS 

It is clear that, with the appearance of km-scale detectors, neutrino astronomy will soon become 
an interesting aid in the study of microquasars and other stellar objects. Here we have focused 
on showing that the current state of the art in neutrino observations can also be used to establish 
useful constraints on some of the hadronic models presented in the literature. 
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Neutrino energy [GeV] 



Fig. 1 Fluxes of neutrinos produced in the framework of different models of hadronic 
microquasars as compared with the atmospheric neutrino background and the upper 
limits for typical sources obtained by AMANDA-II. 
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